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Constitutive equations are proposed for the relaxation of a glassy material in contact with a ther- 
mal reservoir. The dynamics of a single state variable, identified as an effective temperature accounts 
for the emergence of glassy behavior at low bath temperature. Memory effect in cooling/reheating 
experiments are discussed along with several aspects of the response to shear: aging of the response 
to step strains, response to oscillatory forcing, and steady rheology. 
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The purpose of the present work is to introduce a lim- 
ited set of kinetic equations which describe the out-of- 
equilibrium relaxation of a structural glass and its re- 
sponse to shear deformation. It was originally moti- 
vated by recent theories for the plasticity of amorphous 
solids, 1, 2] in an attempt to incorporate glassy relax- 
ation at an elementary level. 0, A quite simple pic- 
ture emerges, which accounts for important properties of 



glassy materials, while its 
classes of complex fluids; 
structural relaxation 




remises may hold for general 
it echoes early theories of 
in a much more system- 
atic framework, and provides theoretical grounding for 
phenomenological rate-and-state equations. [Io| 

We consider, in this study, a piece of material at pres- 
sure P, in contact with a thermal reservoir at temper- 
ature T. We sort the material into subsystems p| (or 
zones) of, say, z molecules, and characterize structural 
disorder by the distribution p of volumes, Vi G [vq, oo], or 
enthalpies (hi = Pvi) of these molecular configurations. 
(i>o is a lower bound imposed by excluded volume ef- 
fects.) At low T, the distribution p evolves slowly, as op- 
posed to the vibrations of the molecular structure which 
quickly equilibrate with the thermal bath. The evolu- 
tion of p is expected to result from transitions between 
metastable basins in phase space, and these transitions 
are supposed to be triggered by conformational changes 
of molecular subsystems (rearrangements). The typical 
size of a rearranging zones fixes z, which is here supposed 
to be independent of T and P. At all times, molecu- 
lar subsystems are in mechanical contact with neighbor- 
ing subsystems, hence likely to exchange energy among 
themselves: very small displacements of the average po- 
sitions of the molecules are sufficient to ensure these 
volume (or enthalpy) transfers. Mediated by acoustic 
modes (phonons), these transfers are thus expected to 
be fast, collective, processes, in contrast with rearrange- 
ments which, at low T, are expected to be rare, space- 
localized events. 

This observation in the basis of the present work: it 
permits to separate these two thermodynamic processes 
(equilibration with the bath v.s. transfers between sub- 
systems) and assume that spatial enthalpy fluctuations 
equilibrate on timcscales at which the system is not in 
contact with the thermal reservoir. At all time, approx- 



imation is made that spatial fluctuations are in "adia- 
batic equilibrium", characterized by a (P,i/)-ensemble: 
the distribution p is determined at every single mo- 
ment by statistical inference: it maximizes entropy un- 
der the single constraint that H = J dh p(h) h is known, 
whence p(h) = -g cxp(— Xh), with the partition func- 
tion Z = J^ =Pvo dhexp(-Xh) = (1/A) exp(— Xho), and 
where A is the Lagrange multiplier that enforces adia- 
baticity. = X/(kX) is an effective temperature, but dif- 
fers from the temperature of the thermal reservoir, unless 
the system has reached thermal equilibrium. This effec- 
tive temperature is a generalization of Edwards': I 
made no assumption here that the system is exploring 
all configurations of some ergodic component at any rel- 
evant timescale but, rather, that, by being very large, 
it samples, at any single time, all possible configura- 
tions allowed by macroscopic constraints. At all time, 

ke(t) = H(t)-Pv . 

The previous Ansatz for p serves here as a basis to 
investigate the relaxation of the system towards equilib- 
rium. At any time, H (or O) suffice to characterize the 
internal state of the material. |16l | Their dynamics is spec- 
ified by energy conservation: it results from the balance 
between thermal exchanges with the reservoir (heat) and 
the work of external forces. Estimating the rate of ther- 
mal exchanges with the bath is a difficult task; it requires, 
in principle, to determine all probabilities p(h — > h') for 
the reconformations of molecular subsystems, leading to: 



^ = / dhdh'p(h -> h') p(h) (h 1 - h) + W e 



(1) 



I argue that a generic form of this equation arises from 
the study of dominant contributions to the integral term. 
Then, I show that the resulting equation accounts for the 
emergence of slow relaxation and glassy behavior at low 
T, with a marked cooling rate dependence of the glass 
transition point. A study of the response to shear in 
various experimental set-up completes my presentation. 

Rearrangements correspond to elementary contacts 
with the thermal reservoir: a fluctuation dh of enthalpy 
requires a transfer of heat Sq = —Sh from the ther- 
mal reservoir, hence occurs with probability exp(5S/k) = 
exp(—5h/(kT)). These thermally activated processes are 
very sensitive to the external constraints imposed on a 



molecular subsystem by its surroundings: for low v, the 
reconformation of the molecules may require to deform 
their "cage" and pay a consequent price in elastic en- 
ergy; for large v, molecules may move freely. Suppose, 
for the sake of simplicity, that the transitions h — ► h! 
are controlled by some energy barrier hb which is a sim- 
ple function of the volume v. Obviously, hb(v) must be 
non-increasing, hence there exist an activation volume 
v a , such that hb(v a ) = Pv a = h a ; beyond this point 
(v > v a ) the notion of barrier breaks down. Two dif- 
ferent families of reactional pathways must be identified: 
when h,h' < h , transitions are activated, and transi- 
tion state theory indicates that the rates must be writ- 
ten: p(h, — > hi) = v exp(— (ft,;, — h)/(kT)); otherwise, 
there is no barrier per se and Monte-Carlo weights read: 
p(h -> ti) = i/nun(l,exp(-(7i' - h)/(kT))). These tran- 
sition probabilities verify detailled balance and are con- 
tinuous functions of either h or hi when they cross h a . 
v is an update frequency, the temperature dependence 
of which might be neglected as it brings only minimal 
corrections. The calculation of the rhs of equation JJJ 
is further simplified by allowing only transitions such 
that h — hi = ±Sq; the integral is separated into the 
sum Ii over h,h! € [ho,h a ], and I2 over the comple- 
mentary domain of integration. It appears that I\ (resp. 
h) is of order O((0 - T) 2 ) (resp. 0(9 - T)) close to 
equilibrium, and is proportional to a factor of the form 
exp(-A/(fcT)) (resp. exp(-A'/(/c9))) when 9 >> T. 
In both these limits, I± is dominated by I2: the dynam- 
ics of 9 is controlled by the transition pathways (found 
with a frequency cxp(— P(v a — Vo)/(k@))) which do not 
involve the crossing of energy barriers. It is therefore suf- 
ficient to restrict our discussion to the situation in which 
Ii — » 0. To fix ideas, I assume that h {v) — * 00 below 
v a and drops sharply at that point. (Such steep decay 
of hb is expected, for example, in the case of hard-sphere 
materials.) 

After integration one gets from equation 0J: 



k 9 = Ei exp 



Ah 



fc9 kT 



exp I — -•—)"• 1 I +<7 T > 

(2) 

with Ei = v Sq and Ah — P(v a — vo). I anticipate here on 
the second part of this paper where the deformation of a 
material in a pure shear geometry will be considered: a 
is the (deviatoric) shear stress, 7 the plastic strain rate, 
whence W cxt = <ry. The results presented here do not 
depend qualitatively on further dependencies of Ei on T 
or 9, which may, however, bring logarithmic corrections. 

In the absence of external stress (a = 0), equation 
admits thermodynamic equilibrium, 9 = T, as its sin- 
gle fixed point. The timescale of the relaxation towards 
equilibrium increases with decreasing T, and at low T the 
system may remain out-of-equilibrium on any observable 
timescale by the mere fact that 9 is small (or tries to). 




FIG. 1: Energy (left) and calorific capacity (right) as a func- 
tion of temperature during cooling and reheating protocol, 
between fcT max = 0.5 and kT m i n = 0.001, for cooling rates, 



10" 3 /fc (solid line), 
line) . 



HT 4 /* (dashed line), 10" 5 /fc (dotted 



When 9 >> T, equation is dominated by: 



k 9 = -Ei exp 



Ah 
fc9 



+ o"7 



(3) 



In this limit, the long-time relaxation of the variable 
9 is easily obtained by integrating JSJl (with 717 = 0): 



f@ e x P [A/i/fe 0] fcd0 = —Eit, and keeping only the 
dominant contribution: k<d(t) 



Ah 



IogCEit/Afc)' FOT fi " 

nite T, this process dominates so long as t « tt — 
Ah exp [Ah/{kT)] /Eh.: t t diverges in the limit T -> 0. 

An important property of the glass transition is that 
the observed transition point depends on the cooling 
rate. [f| Consider the following experiment: enthalpy 
is monitored during a steady cooling from an initial tem- 
perature T max ; after reaching T m - ln , the temperature is 
increased at the same rate back to r max - This procedure 
is reproduced here by integrating equation (J2J for differ- 
ent cooling/reheating rates. The resulting dynamics are 
portrayed figure 0) (left), where energy is shown as a 
function of T . Slowing down is seen as energy does not 
reach equilibrium at low temperature; upon reheating, 
energy does not go back along the same curve. As in the 
experiments, heat capacity is then extracted from this 
data according to the formula, Ct = dH/dT; the result 
is displayed figure (right). The large increase of Ct 
observed upon reheating marks a fictive glass transition 
point, which clearly depends on the cooling rate. 

The rest of this paper is devoted to the study of a ma- 
terial under shear. Like thermal equilibration, the defor- 
mation of a material under an applied shear stress is ex- 
pected to proceed by space-localized rearrangements. 
However, the transition pathways that contributes re- 
spectivelly to elementary shear and elementary thermal 



equilibration processes, are not necessarily identical at 
the molecular level. Therefore they may in general, 
display different activation enthalpies. Without more 
knowledge, we can expect the response to small shear 
stresses to be of the form: 

El 



7 = Eq exp 



Aft' 

'fee 



(4) 



The update frequency Eq normalizes the frequency of ele- 
mentary vibrations, and incorporates corrections specific 
to this mode of deformation. Any dependency of Eq on 
or T is neglected. (Eq/Ei and Ah/ Ah' are expected 
to be of order 1.) 

With equations @ and (0) in hands I now study 
three experimental protocols: aging in the stress-linear 
response to step strains, response to oscillatory forcing, 
and rheology of steady deformation. 

Aging is characterized by monitoring the response of a 
material to a perturbation at various times after an initial 
quench from a high temperature equilibrium state. |flll2^ 
The procedure is as follows: At time t = 0, the system is 
suddenly quenched from an equilibrium state (6o = To) 
to low temperature. Relaxational dynamics take place 
from then on, and we study the situation when >> 
T: equation J3Jl leads at long time to a time-logarithmic 
decay of enthalpy. At time t w , a small strain Se is applied; 
it provokes a quasi-instantaneous elastic response a = 
2/j.Se, followed by a slow occurence of plastic deformations 
and stress release as: a = — For small Se, the term 
7<7 in equation Q is negligible at all times. 

At short times, stress decays exponentially, with a 
timescale, tq — (l/(2fiE )) exp [Ah' / (k<d w )], which de- 
pends only on Q w = Q(t w ). At long-times, stress release 
is controlled by the slow relaxation of 0, and integrates 
as: 



G{t,t v 



*(*) 

a(t w ) 



exp [A 



with = 1 - % and A = Zfi. (f^ . If 
Ah > Ah', (3 > and stress undergoes KWW relax- 
ation. The crossover between the short-time exponential 
relaxation and the long-time streched exponential is il- 
lustrated figure A plateau appears in the crossover 
region: this typical pattern emerges here solely from a 
change in the dynamical regime of 0. Relaxation curves 
for different ages t w are shown figure 01 (top). If fit- 
ted by exp(— t /r(t w )), it is easy to check that the ap- 
parent relaxation timescale r(t w ) grows like, t ^, with 
a = 1 — j3; this phenomenon, usually refered as sub- 
aging is here a mere artifact of the lin-log representation: 
this type of scaling - which are very often used to treat 
experimental data - may be very misleading and hide, 
as it does here, the true aging behavior. If Ah = Ah', 
[3 = 0, wc find G(t) ~ 1/t. If Ah < Ah', (3 < and 
stress assumes a non-vanishing asymptotic value after 




logioO) 



FIG. 2: Parameters are E = fj, = 1, Ah = 1.5, Ah' = 1, 
0(0) = 10, and t w = 0. Top: I/O as a function of log 10 (t) for 
E-t = 1,20,40,80,100 from right to left. Bottom: relaxation 
spectra with E\ increasing from left to right. 



relaxation, which increases with the age of the sample: 



G(t — > oo ) ~ exp 



\A\t\ 



Another important experimental protocol consist in 
measuring the in-phase and out-of-phase moduli, G' and 
G" , when the material is forced at time t w by an oscil- 
latory shear, j(t) — 70 sm(oj (t — t w )). In practice, the 
measurement is performed by integrating the response a 
fixed number of periods n, for each frequency to. Since the 
material ages, every point on the spectrum depends on 
the whole history of the sample prior to its measurement. 
The resulting G' and G" are shown figure (bottom) . 
They clearly show an a relaxation peak which, here is 
solely the signature of the slow evolution of 0. This peak 
is obtained for large intervals of the parameters, under 
the conditions that /1 is somewhat large (which, indeed, 
is reasonable, given the usual values of elastic constants); 
it weakens with the increasing age of the material. 

To conclude this study, let me consider the regime of 
stationary deformation of the material, at temperature 
T, under a constant shear rate 7. The steady state value 
of the variable is the solution of: 



7 = Eq Ei exp 



Ah + Ah' 



Sq Sq 



Ah 

Sq 



For Ah — Ah', one finds, 777 
For Ah j£ Ah', and fc0 >> Aft + Aft', this solu- 
tion is dominated by the contribution, = kT/(l 4 



4 




FIG. 3: Top: Relaxation modulus for Eo = Ex = 
H = 1, Ah = 1.5, Ah' = 1, 6(0) = 10, and t w = 
10, 10 2 , 10 3 , 10 4 , 10 5 , 10 6 , 10 7 . Bottom: Normalized frequency 
spectra for fi = 100, E = Ex = 0.01, Ah = 1.5, Ah' = 1, 
t w = 10, 6(0) = 10. 

(kT/Sq) log(l -^/(E x Eq))) which, for small 7, indi- 
cates that the system behaves as a newtonian liquid with 
a viscosity, rj = exp [Ah'/kT] /E Q . For kT ^ fc9 << 
A/i + A/i', the solution is dominated by its T — > value, 
which leads to a power relation between stress and strain 
rate: a ~ 7", with n = ■ 

This study suggest that, in the low temperature 
regime, there exists a direct relation between the expo- 
nent (3 (and a) of the KWW relaxation, and the power 
law rheology displayed in steady shear. It establishes 
a relation between two highly dissimilar experimental 
situations - linear response in the aging regime, and 
steady shear deformation - in direct relation to the ra- 
tio k = Ah J Ah! between activation barriers associated 
respectivelly to shear and energy relaxation processes. 
k is expected to depend solely on geometrical aspect of 
the material: shape of the molecules or polydispersity of 
colloidal particles. The existence of such a relation can 
thus be tested experimentally by considering families of 
closelly related materials in order to allow some variation 
of the exponents. 

The equations proposed here are exaggeratelly sim- 
ple. Their interest lies in the possibility to account for 



a whole set of properties commonly associated with the 
glass transition, in a over-simplified framework, in com- 
parison to other studies of rheology. 0] Various aspects 
of these questions certainly deserve further studies. In 
particular, a derivation of the transformation rates would 
bring out much needed information. One may, however, 
hope that out-of-equilibrium thermodynamics can be de- 
fined for structural glasses in terms of a few rate equa- 
tions for the appropriate set of thermodynamic quanti- 
ties. 
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